Letter 



Typeset with jpsj2.cls <ver.l.l.l> 



Upper Critical Field of the 3 Kelvin Phase in Sr 2 Ru04 
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The inhomogeneous 3 Kelvin phase is most likely a superconducting state nucleating at the 
interface between micrometer-sized Ru-metal inclusions and Sr2Ru04 above the bulk onset of 
superconductivity. This filamentary superconducting state yields a characteristic temperature 
dependence of the upper critical field which is sublinear, i.e., H C 2(T) oc (T* — T) 7 with 0.5 < 
7 < 1 (T*: nucleation temperature). The Ginzburg-Landau theory is used to analyze the 
behavior of the nucleated spin-triplet phase in a field and the characteristic features of H C 2 
observed in the experiment are explained based on a two-component order parameter in the 
presence of a filament of enhanced superconductivity with a finite width. 
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Sr 2 Ru04 plays an exemplary role among unconven- 
tional superconductors as a realization of spin triplet 
pairing in a quasi-two-dimensional (2D) Fermi liquid, 
with some similarities to superfiuid 3 He. 1-3 ) Experi- 
ments provide strong evidence of a superconducting state 
with in-plane equal-spin pairing 4 ) and violation of time 
reversal symmetry. 5 ' This uniquely identifies the pair- 
ing symmetry to be that of a chiral p-wave state, analo- 
gous to the A-phase of 3 He: d(k) — An£(fc x ± ik y ). 6 ' 7) 
This is a lucky case in many respects. We mention only 
a few points here. (1) Broken time reversal symme- 
try is responsible for unusual magnetic properties. (2) 
The order parameter consists of two complex components 
t] = iWx-> Vy)-> the only case among all possible triplet pair- 
ing states, for tetragonal crystal symmetry, which cor- 
respond otherwise to one-component order parameters. 
Thus, we may write 



d(k) = z(r) x 



Tjyky 



1 



(3) It gives rise to unusual vortex physics, including a 
square vortex lattice and anomalous low-temperature 
flux dynamics. 8 ' 9 ' (4) Chiral gapless subgap quasipar- 
ticle states appear at the surface. 10 ' n ' 

For these properties the important feature is the de- 
generacy of the two order parameter components, which 
is guaranteed by the tetragonal symmetry. It was sug- 
gested that this degeneracy lifted by symmetry lowering 
would lead to two consecutive phase transitions similar to 
those in the heavy Fermion superconductor UPt3. One 
of the way to realize this is the confinement of the su- 
perconductor in a narrow filament which has by geome- 
try a symmetry lower than tetragonal. 12 ' Such kind of 
filamentary superconductivity is most likely realized in 
an inhomogeneous superconducting phase dubbed the " 3 
Kelvin phase" of Sr2Ru04. 13, 14 ' This phase appears at 
nearly double the bulk transition temperature T c = 1.5 
K and is associated with the presence of micrometer- 
sized Ru-metal inclusions in the otherwise very clean 
Sr2Ru04. It has been suggested that this phase nucle- 
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ates at the interface between the Ru-metal and Sr 2 RuC>4, 
where the critical temperature is larger possibly due to a 
locally enhanced density of states and modified electron- 
electron interactions. 15 ' In such a case the supercon- 
ducting state appears at a temperature T* > T c in a re- 
stricted region of lower symmetry. This superconducting 
state has a single order parameter component and does 
not violate time reversal symmetry. The component cor- 
responds to the p-wave superconducting state with mo- 
mentum along the interface, i.e., n • rj = where n is the 
interface normal vector. This filamentary phase yields 
several unusual properties. Unlike in a conventional s- 
wave superconductor, the transition to the bulk phase 
is not merely a matter of percolation, but represents a 
real (time reversal) symmetry-breaking transition. This 
corresponds to an additional second-order phase transi- 
tion. 16 ' Moreover, this system may constitute a complex 
intrinsically phase frustrated superconducting network. 

An important way of probing the filamentary phase is 
the observation of the nucleation in a magnetic field, i.e., 
the upper critical field H c i- The confinement to a nar- 
row filament yields a sublinear temperature dependence 
H c2 {T) cx (T* - Tp where 7 = 0.5 is close to T* in 
contrast to the linear behavior for the bulk H C 2. 15, 17, 18 ' 
In view of experiments showing exponent 7 lying be- 
tween 0.5 and 1 in qualitative agreement with the ex- 
pectations, 19 ' we would like to reanalyze the behavior of 
H C 2 in the filamentary phase. 

Our analysis is based on the Ginzburg-Landau model 
of an infinite planar interface, as introduced in ref. 15. 
It was suggested that the locally enhanced T c at the in- 
terface is the result of a local lattice distortion mainly by 
RuOg-octahedra rotation around the z-axis. This gives 
rise to reduced hopping matrix elements such that the 
Fermi velocities decrease, increasing the density of states 
with an additional (Stoner) enhancement of the uniform 
spin fluctuations. The extension s of the distortion is 
of the order 100 A and, 20 ' thus, is much shorter than 
the coherence length £. On the other hand, the size of 
the Ru inclusions is ~ 1 — 10 /im, large compared to 
the coherence length. 13,14 ' Hence we consider here an 
infinitely extended interface separating two half spaces. 
The purpose of the present study is to analyze the nucle- 
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ation of the filamentary unconventional superconducting 
state in a field. Therefore, we will further simplify our 
model. The system is taken to be symmetric at the in- 
terface so that the material on both sides is identical. 
This deviation from reality has only a minor influence 
on the qualitative behavior of H c2 . Because the specific 
orientation of the interface is not so essential, we choose 
a concrete example in which the normal vector lies in the 
a — b plane and points along the x-axis. The interface 
location is then defined by x — 0. Our model is given 
by the standard Ginzburg-Landau free-energy functional 
for the above p-wave order parameter (i] x ,ri y ). The lo- 
cal enhancement of superconductivity at the interface is 
introduced by an additional interface term: 15 ) 

F = [ dx{a(\r, x \ 2 + \ Vy \ 2 )-^6(x)(\ Vx \ 2 + \ Vy \ 2 ) 

+ ^i(|A^| 2 + |A,%| 2 ) + K 2 (\D yVx \ 2 + \D X% \ 2 ) 
+ [K 3 (D x r) x )*(D y r]y) + K 4 (D y r) x )* (D x r) y ) + c.c] 

+ K 5 (\D z r) x \ 2 + \D z r) y \ 2 ) + ^-(V x A) 2 j. (2) 

This term is localized to the interface by a delta func- 
tion, since the extension s is small. Although we will 
see later that for a higher field this extension can be im- 
portant, we will ignore it for the moment. Moreover, 
D = V + i(2e/hc)A and a = (T-T c )/T c , where A is the 
vector potential and T c = 1.5 K for the bulk Sr 2 Ru04. 
Ki (i = 1,2,3,4 and 5) are coefficients for the gradient 
terms ( Ki/3 = K 2 = K$ = K4 for a cylindrical Fermi 
surface) . We define £ = y/Ki as the characteristic length 
scale of the order parameter in the bulk material. For 
the discussion of the nucleation of the superconducting 
phase, we can neglect the quartic terms in the Ginzburg- 
Landau theory. 

The variation of F with respect to r\ x and r\ y in a zero- 
field, leads to two decoupled differential equations in the 
form of Schrodingcr equations: 



[-K 1(2) d 2 x - t]crS(x)] r] x{y) = -an. 



•x( v )- 



(3) 



The "lowest energy" solution corresponds to a bound 
state for an attractive delta potential (locally enhanced 
superconductivity) leading to an eigenvalue a which de- 
termines the nucleation temperature. The order param- 
eter components shows an exponential decay towards the 
bulk region: 

Vx(y) = exp(-|z|/6r(2/)), with £ x(j/) = 2K 1(2) /£<7 (4) 

and the transition temperatures for r\ x and r\ y are ob- 
tained from the eigenvalues a = (T* — T c )/T c > 0: 



and T;=T c (l + ^g). (5) 



Both are enhanced by the a term and T* > T*, since 
K\ > K 2 . Hence the first nucleation of superconductiv- 
ity occurs in the ^-component whose nodes lie perpen- 
dicular to n. 

We study first the case of fields in the x — y plane. 
The upper critical field H c2 for the nucleation of super- 
conductivity is higher for H _L n (|| y) than for H \\ n. 
For H || y, the components r\ x and r\ y remain decoupled 



at the nucleation point. H c2 is determined by the in- 
stability of rj y . We use A = (0, 0, — Hx) leading to the 
equation 



[-K 2 d 2 x - &5{x) + 



4K 5 x 2 



\Vv 



-ar) y 



(6) 



for rj y , where Ih = \Jch/(eH) is the magnetic length. 
For low fields (long Ih), the harmonic potential term is 
a weak perturbation to (3), so that i] y = cxp(— |x|/£ y ) 
remains a good approximation. Substituting it into (2) 
and integrating over x we obtain the second-order term 



(7) 



H 



whose zero determines the instability. Thus F = yields 
H c2 : 



H C 2 = — 1 



K 5 e 



IT* — T 



(8) 



as found previously. 15, 17, 18 - ) With increasing field, how- 
ever, Ih becomes shorter and the harmonic potential 
term in (6) becomes a larger correction. Thus the exten- 
sion of i] y along n shrinks due to the decreasing cyclotron 
radius. A good approximation to the ground state of the 
"Schrodinger equation" (6) is obtained by a variational 
ansatz, which captures the basic behavior in a simple 
way: 



r] y = exp( 



\x\/Z v )eM-VK 5 /K 2 x 2 /l 2 H ). (9) 



Here exp(— y/K 5 /K 2 x 2 /l%) describes the asymptotic be- 
havior for distances far from the interface, while the ex- 
ponential part gives the correct boundary conditions at 
the interface. Again we substitute i] y into (2) and in- 
tegrate over x, so that setting F = leads to H c2 (T) 
(Fig. 1). We may approximate the low-field range, by a 
power law H c2 cx (T* — T) 7 . We find the best approx- 
imation to the variational solution for 7 = 0.62 in the 
range 0.9T* < T < T* = 2.8 K. This value compares 
well with recent experimental findings of 7 = 0.7 — 0.75, 
which is indeed sublinear. 19 ) The exact square-root be- 
havior in the limit of very small fields is experimentally 
difficult to observe. The limitation of this behavior is 
given by £ -C Ih, i-c, H <C H* with a characteristic field 
H* = 2^/K 2 /K b (<S> /2ir)l£ 2 - IT. Moreover, in Fig. 
1 we observe a deviation from our variational solution 
for T < 2K. This is partially due to the limited validity 
of the Ginzburg-Landau theory which only extends to 
the region close to T* . Furthermore, the suppression of 
H c2 is likely related to a limiting effect (analogous to the 
paramagnetic limiting) which has also been observed in 
the bulk H c2 for fields in the basal plane. The discussion 
of this high-field behavior lies beyond our model and our 
scope. 

For H || n, we use A = (0,0, Hy) leading to the 
Ginzburg-Landau equation for rj y : 

[-K 2 d 2 x K x d 2 y + &5{x) + ^^]vy - -any. (10) 
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Fig. 1. Temperature dependence of the upper critical field for an 
in-plane field. The solid line denotes H \\ y. For a low magnetic 
field, H c 2 has a square-root dependence H c 2 oc (T* — T) 5 , 
which is plotted as a dashed line for H \\ y. The long-dash 
linear line denotes H \\ x. We choose the following parameters: 
KrJKr = 1/500 and K 1 /(cH/2e) = ^ 2 /(ch/2e) = 20. cr is given 
to fix the transition temperature T* = 2.8 K. Circles (stars) 
denote the experimental data for down (up) sweep of field or 
temperature, 19 ^ showing a hysteresis in H c 2 below 1 K. 



Since the x- and y-dependences factorize, we obtain the 
relevant solution: 



Vy 



eM-\x\/^y)eM-VK^/Kiy 2 /l 2 H) , (") 



and H c2 oc T* — T as shown in Fig. 1 (long-dash line). 
The nucleation field for this direction is obviously lower. 
It is clear that the observed H c2 is due to interfaces that 
lie parallel to the applied field (H _L n), which yields 
the highest nucleation field. 

Now we turn to fields parallel to the z-axis assuming 
simultaneously H _L n. In this case the two order pa- 
rameter components are coupled. We choose the vector 
potential as A = (0, Hx, 0) so that the free energy is 
expressed as 



F 

fx(y) 



/GO 
dx(f X + fy + f xy ), 
-oc 



= a\Vx(y)\ + Kl(2)\dxVx(y)\ ~ fr S (. x )\Vx(y)\ 



+ - 



4K 



2(1)3! 



I 4 



\Vx(y)\ 



2(K 3 + Ki) 



I 2 

1 H 



{ir) x d x rf + irj y d x r]* + c.c.}. 



(12) 

As in the case of H || y, we introduce a variational form 
for the order parameters: 



V v (x) = c y{x) exp(- 



^)eM-jK m /K 2il)JT ), 

?j/(x) v l H 

(13) 

where C x r y \ are coefficients to be determined in order to 
maximize the nucleation temperature for a given mag- 
netic field. We can integrate the free energy analytically 
and determine the upper critical field explicitly We plot 
the result in Fig. 2. Our result well reproduces the ex- 



Fig. 2. Temperature dependence of the upper critical field for 
H || 2. For quite low magnetic fields, H c 2 has a square-root 
dependence, which is plotted as a dashed line. The parame- 
ters are the same as in Fig. 1. Circles denote the experimental 
data. 19 ) 



periment above 0.9T*. In this region, we find a exponent 
7 = 0.76 fit to a power law, which agrees well with the ex- 
perimental result (7 = 0.72). 19 ) In the vicinity of T* H C 2 
exhibits again a very-low-field square-root dependence as 
in the case of if || y, which is plotted in Fig. 2 as a 
dashed line. In this case, H c2 deviates even more rapidly 
from the square-root behavior as temperature decreases 
than the in-plane fields. The reason lies in the lower 
characteristic field H* = 2^ ' K 2 j ' K 1 {^ /2tt)( 1 - 0.05T. 

In Fig. 2 the experimental data show a pronounced up- 
turn for T < 2K, opposite to the trend for in-plane fields. 
We would now like to discuss the origin of this behavior 
by extending our model. So far we have assumed that the 
extension of the region with enhanced superconductivity 
is negligible and is well described by a delta function. 
We replace, however, now the delta function in eq.(2) by 



6(x) 



1 



y/TTS 



cxp[— {x/ s)" 



(14) 



where s represents the width of the interface region 
of enhanced transition temperature. In the introduc- 
tion we speculated that this region is characterized by 
an increased density of states or, equivalently, by a 
decreased Fermi velocity. Since the coefficient of Ki 
(i = 1,2,3 and 4) is connected with the Fermi veloc- 
ity in the x — y plane (fQ cx v 2 F /T 2 ), we introduce in 
addition, the following spatial dependence in Ki. 

Ki -» Ki{\- Aexp[-(x/s) 2 ]}. (i = 1, 2, 3 and 4) (15) 

with A (0 < A < 1) as a parameter. Assuming the same 
variational order parameter form we integrate the free 
energy analytically. The resulting H c2 is shown in Fig. 
3. In particular, we observe the onset of an upturn of 
H c2 for H || z at approximately 2 K, which compares 
well with the experiment. 

This feature originates from the fact that with shrink- 
ing magnetic length the region of nucleation becomes 
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Fig. 3. Temperature dependence of the upper critical field for H \\ 
y and z. Solid lines are the results for a finite s. Dashed lines 
represent the result for s = 0, which are same as in Figs. 1 and 
2. We choose the following parameters for a finite s: s = 0.3£, 
A = 0.7, K 5 /K 1 = 1/700 and K 1 /(ch/2e) = £ 2 /(ch/2e) = 24. a 
is fixed to give the transition temperature T* = 2.8 K. We used 
the same experimental data as in Figs. 1 and 2. 19 ) 



increasingly confined into a narrow layer, where we find 
an enhanced transition temperature as well as a shorter 
local coherence length. Both act to increase the criti- 
cal field. Thus the onset of the upturn is expected when 
Ih ~ s. The fit to the data with our variational approach 
yields s ss 200 A. Note that this kind of upturn behavior 
is not expected for in-plane fields, since K 5 , that deter- 
mines the coupling, would not have the form (15) as it 
depends on the z-axis Fermi velocity that would not be 
significantly affected by RuOe-octahedra rotations. 

In summary, the discussion of the upper critical field 
for the filamentary superconducting phase exhibits sev- 
eral length scales to be taken into account, which are 
the effective magnetic length Ih, the coherence length £ 
and the extension s of the interface regions. If the ef- 
fective magnetic length Ih is much larger than £ and s, 
we find that i/ c2 would follow the square-root power law 
(T* — T) 1 / 2 . Once the field starts to shrink the extension 
of the nucleated order parameter, we encounter a devia- 
tion from this behavior and a power law fit would yield 
a different power law. Our discussion shows that in the 
case of the 3-Kelvin phase the square-root behavior oc- 
curs in a very limited range of low fields which is hard to 
analyze. Finally, if the field becomes sufficiently strong 
to confine the nucleating order parameter in the inter- 
face region s a relative increase of the upper critical field 
is possible. In the 3-Kelvin phase this is observed for 
the field along the z-axis. However, it has to be noticed 
that an additional important enhancement factor is the 
coupling of the two order parameter components. This 
latter effect is due to the Zeeman coupling of the mag- 
netic field to the orbital magnetic moment of the Cooper 
pair for the chiral p-wave phase. 

The comparison of our theory with the experiment 
shows that we are in principle able to fit the experimen- 



tal data. However, this aspect has to be viewed with 
care, as the Ginzburg-Landau free energy is expanded 
at a temperature around T* and has therefore limited 
quantitative reliability. Moreover the model was in many 
respects simplified. Nevertheless, the ability to repro- 
duce the qualitative features are convincing, we believe. 
One obvious problem is the limiting behavior for in plane 
fields which seems to be present in both the bulk and 3 
Kelvin phase of Sr 2 Ru04. 

The upper critical field may be considered a strong 
confirmation of the assumption that the 3 K phase 
Sr 2 Ru04 is due to the local enhancement of the su- 
perconducting transition temperature at the interface 
of Ru-inclusions and Sr 2 Ru04 . Thus, we expect that 
Sr 2 Ru04 has two consecutive phase transitions due to 
the symmetry lowering by the Ru inclusions. This fact 
still remains to be experimentally verified. It could then 
be added to the other convincing evidence for chiral p- 
wave pairing in Sr 2 Ru04. 
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and the MaNEP project of the Swiss National Fund. 
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